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CN ■ Abstract 

>. 
^D . Propagation of gluons in the confining vacuum is studied in tlie 

framework of the background perturbation theory, where nonpertur- 
\^ ' bative background contains confining correlators. 

Cp . Two settings of the problem are considered. In the first the con- 

fined gluon is evolving in time together with static quark and anti- 
quark forming the one-gluon static hybrid. The hybrid spectrum is 
""Oj' calculated in terms of string tension and is in agreement with earlier 

JL . analytic and lattice calculations. In the second setting the confined 

Q I gluon is exchanged between quarks and the gluon Green's function 

>-^ ■ is calculated, giving rise to the Coulomb potential modified at large 

distances. The resulting screening radius of 0.5 fm presents a seri- 
k>( I ous problem when confronting with lattice and experimental data. A 

;h ' possible solution of this discrepancy is discussed. 

1 Introduction 

Gluons are known to be confined but this property is never taken into account 
in the Standard Perturbation Theory (SPT). As an argument one refers to the 
small distance (high momentum) domain where SPT is assumed to be valid. 
Beyond this domain the SPT displays unphysical singularities, and moreover 
the very notion of gluon should be properly defined. This can be done in the 
framework of the Background Perurbation theory (BPT) pp. The formalism 
of this kind where the background is assumed to be nonperturbative with 
confining properties, was developed in j2l E] • 
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One immediate consequence of this new BPT is that all gluons are con- 
fined and moreover the unphysical singularities of SPT (Landau ghost poles 
and branch points as well as IR renormalons) are removed from the theory 

urn- 

The confined gluons can form several types of systems: glueballs j3], 
hybrids E] and gluelumps [7|. The analytic calculations in the quoted 
papers predict spectrum which in all cases is in good agreement with lattice 
data, and has a very simple form depending only on string tension a and as- 

The main subject of the present paper is the study of the gluon exchange 
interaction between Q and Q when gluon is confined, which can be called 
the confined Coulomb interaction V*{R). One can expect that at small 
distances, say R < 1 GeV~^, the confined Coulomb potential coincides with 
the standard Coulomb potential K(-R) = —C2 "^^ , C2 = ^n ■ 

At large R the confined gluon is expected to produce the screening of the 
Coulomb interaction. At the first glance the screening mass should coincide 
with the lowest hybrid mass of the static hybrid. However as will be shown 
below, this naive expectation fails, since gluon propagation between quarks 
goes not only in time (where hybrid mass is in proper place) but also in 
distance R (where in addition asymptotics of wave functions enters). As 
a result one obtains a more complicated behaviour which we display both 
numerically and analytically. 

The plan of the paper is as follows. In section 2 the general formulas 
of BPT are written down and approximations are discussed. In section 3 a 
simple toy model is suggested to illustrate method and possible qualitative 
outcome. 

In section 4 the method is applied to calculate the static hybrid Green 
function and spectrum in a way different from used before, in |E|. In section 
5 results of previous section are used to calculate the Green function of the 
confined gluon exchanged between static quarks, and the resulting screened 
Coulomb potential. Physical consequences and prospectives are discussed in 
the concluding section. 



2 Dynamics of a confined gluon. Generali- 
ties. 

In the Field Correlaror Method (FCM) the dynamical picture of a confined 
gluon is simple and self consistent: the gluon (its corresponding field is a^) 
moves in the strong and disordered vacuum field B^ characterized by the 
correkators of the field F^y{B)^ so that the total gluonic field A^ can be 
written as 

A^ = B^ + a^. (1) 

Here the problem of separation of A^ into B^ and a^, and that of double 
counting is resolved technically by the use of the so-called 'tHooft identity 
[3J and one can integrate and average over both Z^-B^ and Dttfj,, so that the 
partition function is 

Z = ^JdB^Z{J,B), Z{J,B)= J Da^exp{-S{a + B) + J{a + B)) (2) 

and 5* is the standard QCD Euclidean action including ghost and gauge-fixing 
terms (see |21IS] for details). 

In what follows we shall be interested in the gluon propagation in the 
field of static (confined) quark Q and antiquark Q. The starting point for 
the gauge-invariant study of this process is the total Green's function of QQ 
system, which is proportional to the Wilson loop: 

{W{A)) = {{W{B + a)))B,a = m{B + a))a)B, (3) 

WiA) =tr P exp ig J A^dz^. (4) 

The Wilson loop is assumed to be the closed rectangular contour R x Tq 
in the {xi,Xi) plane. 

One may expand in ga^^ keeping B^ intact and this will give the Back- 
ground Perturbation Theory (BPT) series ([1112], see [S] for details), the first 
terms being 

{W{B + a))B,a = {W^''\B))b - g^{W^^\B, x, y)G,,ix, y))Bdx,dy, + ... (5) 

For the chosen contour C and the Feynman gauge of field B^ and FFSR for 

Gjjivi W^'^'^ and G^,^ can be written as 

iy(^) = C2(/)$y(a;,y)<l>W(y,a;) (6) 
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G,,{x,y) = J^^ ds{Dz),ye-'' i^<^^''''^\x,y)Pexp{2g J^ Fp^{z{r))dz 



tz]\ 

) f3a,f3'a' 

(7) 
where $*^^^(x, y) are the future/past pieces of the Wilson loop W^^\B) = 
Pexp{ig § B^dZfj), obtained by cutting it at points x and y. Note that the 
adjoint phase factor $(°"^-'^(x, y) = exp(ig Jy B^dz^) from G^j, which can be 
written as the product 

$('^'^^)(a^,y) = '^/3a(l/,a;)<l>;3w(x,i/) (8) 

the total construction produces two closed Wilson loops (see ref. |2l 13 El for 
pictures and discussion). 

$(-)$M)$(+) = wy\x,y)W^+\y,x). (9) 

Note the subscript cr in (jHl) which implies the color-magnetic spin factor (the 
last factor on the r.h.s. of ((Zj)) entering in all Wilson lines including ^^'^'^^\ 
The averaging of © over fields -B^ can be easily done at large Nc 

{W!,I^G,,)b = G^'\wi-\x,y))B{Wi-^\y,x))B, (10) 

where G^^^ denotes the integral / ds{DZ)xye~^ . In the FCM one obtains for 
(W) the area law behaviour at large distances R,T ':^ \ 

{W) = exp(-a5™„) (11) 

and A is the gluon correlation length j5|, characterizing the fall-off in x of 
the field correlator D{x) ~ {trF{x)^F{0)), and Smin is the area, which is 
assumed to be minimal for the given contour C. 

Note that at large distances the spin factor does not contribute to the 
area law [4-7] and therefore the subscript a in (fTT|) is omitted. 

To describe Smin one can parametrize first the trajectory z^(t) of the 
gluon between the initial point x and the final point y (both on the contour 
C). 

z^:{zi{t)=^{t), Z2{t), Z3{t),z^{t)=t), h\t) = zl{t) + zlit) (12) 

We choose the Nambu-Goto ansatz for the minimal area surface (or rather 
for its increase over the plane area 5*0 = RT) 

AS = ASi + AS2, ASi = r dt f df5J{w,w[Y - wfw'^ (13) 

Jyi Jo 
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where Wi^^r, P) : Wi4 = t,Wi2 = (1 — /9)(R, 0) + /3z. Note, that in our case 
two different processes can be initiated by the exchanged gluon: 

i) points X and y are at X4 = Tq, x = (^, 0, 0) and 2/4 = 0, y = (-|, 0, 0), 
the situation which is insured on the lattice by the insertion of plaquettes at 
these sides of Wilson loop. In this case one obtains the hybrid excitation of 
the Wilson loop, and this form was used before to compute hybrid spectra 
analytically [HIEI and on the lattice p. 

ii) the points x and y belong to the trajectories of Q and Q respectively, so 
that G^u{x,y) = G44{R,T), T < Tq, describes the propagation of Coulomb 
gluon between the quarks. 

In what follows we shall study both cases using for that the final form 
resulting from (pi. p. (ITT |) . (P|) 

{Gm{x, y))B = G{x, y) = / rfs(Dz),,e-^-'^^^ (14) 

To treat the awkward roots in AS*, Eq. (fT^ . one can introduce the einbein 
parameters z/(t) and z/(t) JU] obtaining in the small oscillation limit 

G{x,y)= ds{Dz),yDuDue-^^-''^^^''^''^ (15) 

where 

(16) 
For X = (f , 0, 0, To) ; 2/ = (f , 0, 0, o) one will have from (HH) the static 
hybrid spectrum, to be compared with previous calculations, and for the case 
ii) one can define the generalized Coulomb interaction as 

fTo fTo pTo Qf. _|_ y fTo/2 

/ dxj dy,G{x,y)= d^^ G{R,T)dT ^T^V^{R). {17) 

Jo Jo Jo 2 J-To/2 

In the limit a ^ one obtains the free gluon propagator G^^\x,y) = 
4^2 (-^_ )2 in case i) and the Coulomb interaction V*{R) = K(-R) = —(^2^ in 
the case ii) . Our purpose in what follows is to obtain modification of these 
results for nonzero a. To grasp the idea we shall start with a toy model as a 
warm-up. 



3 A simple toy model for the confined gluon 

Consider a gluon propagating from static quark Q with coordinate (0, 0) to 
antiquark Q{T, _R, 0, 0). 

The world sheet of QQ system with the string from Q to Q sweeps the 
strip in the (x4,a;i) plane, and one expects that gluons are confined dynami- 
cally to some region around this strip. This means that the running away of 
gluons from the plane (x4, Xi) is damped by some function, which we take in 
the form of the confining "potential" V , 

V{7.) = ^{zl + zl). (18) 

The Green's function of the gluon can be written in the FFSR ^T] 

rOO r-s 

G{T,R)= / ds{Dz,)oR{Dz^)oT{Dz2)oo{Dz,Ue-''-Jo'"'^ (19) 
Jo 

where K = I J^ dr {^)\ and (D^,),, = WLi '^e^^'^'^'^''^ (^^). 

The integration in Dzi factorizes and can be performed immediately, with 
the result (see Appendix 1 for details) 

G{T,R) = -— / 4y,(u;.)exp(-^), ^{t) = -—; x' = R' + T\ (20) 
loTT^ JO s^ As smht 

The integral (|20|) can be estimated by the stationary point method and 
one obtains 

G(T,m^tk^,m4^: '"'] • (21) 

47r^a;^ [ v^e"^, t > 1 

One can see that at large x^, x'^uj ^ 1, there appears in ()21|) the damping 
factor signalizing the mass gap m = y/u , i.e. the confined gluon acts at large 
distances as a massive particle, while at small distances, x^ -^ 0, it behaves 
as the ordinary unconfined gluon. When cu — > one recovers from ()17|) the 
standard Coulomb interaction: 



PCX) 

V,{R) = -g'C^im G{T,R)dT 
Jo 



asC2if) 



R ' 

and in the opposite limit, uR^ ^ 1, Vc{R) is multiplied by the factor 
. ^2.^1/4 exp(— (cj-R^)^/^). This screening factor is equal 1/2 at i? ^ J^ and 
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decays exponentially at large R, as it is commonly expected. In section 5 a 
more complicated behaviour will be obtained in the realistic case when gluon 
is confined by the string world-sheet. To this end we develop in the next 
section the necessary formalism for the hybrid Green function. 

4 Spectrum of the confined gluons. Static 
hybrid 

In this section we shall calculate the gluon Green's function in the hybrid 
situation, i.e. when boundary conditions are given in i) below Eq. (jT^. and 
one can identify T = Tq. One can introduce the einbein variable fi{t) as in 
[TUj (see Appendix 2 for details) and write 



G{x, y)= j ^DvDve-^G^{R, T, u, u, fi) 


(22) 


where we have defined 




Y- ^^'dt\ "" Tf- 1 -)dt 1 jp' r "^^ 


(23) 


^ Jo 2^^ ' 2 Jo ^ ' ^ ' "^^ Jo 2iu + uy 


G, = |(D^.).,e-r(t^?+^^^^^^i)'^*-^r[(--^)^+'^^l^^ 


(24) 


i? = i?^, and we take x=(f, 0,0, O), y = (f , 0, 0,t) , J, = f i/^. 





The integration over (D^z) can be immediately done using the standard 
path integral formula, see Appendix 1. 

_ , fiuji y/2 (^ _^ j^ _^ J2)uj± 



27rsinhci;iT/ 27r sinh(ci;_LT) 

exp I ^f ^^]^^ ' \l\ cosHu;,T) - 1)1 (25) 

^1 2sinhtUiT 2(z/ + z/)2 ^ ^ ^ M ^ ^ 



where ui = ,/-^, u!± 



vv 



The next step is the stationary point analysis of Eqs. (j^ . (^^ with 
respect to variables u, z/, fi. We relegate the details of this analysis to the 
Appendix 2 and here only quote the result. Minimizing the expression in the 



exponent at (|25|) with respect to {u — v) one obtains that v = v ^ and the 
stationary point in z/, i/ = z/q is found from the equation 

One can distinguish two cases: a) aR^ ^ 1; b) aR^ <^ 1. 

a) In the case a) taking into account that /io ^ ^/ct (which will be con- 
firmed afterwards), one has 

R (0) f^^V (0) Vl2 ,_. 

'' = 2 ' "^ ^ [Hr) ' "- ^ ^ (2^) 

and Eq. ()22j) with u = u = uq assumes the form for large uT 



G{x,y) = jDf,exp[-{aR + Jl^ + | + l^ ^)T] (28) 

The stationary point /x = /xq is found from the exponential of (PH|) to be 



"' \R 
and the resulting static hybrid mass at large R is 

3/(T\l/3 ^ 



Mf,y,„d{R) = ^^ + 2 U J ^ ~^ ^^^^ 

The second term on the r.h.s. of ()29|1 is the characteristic " R~^^^ law" 
for large i? hybrid excitations, studied in 0, and one can distinguish the 
longitudinal and transverse branches of spectrum with higher excitations 
generated by sinhcuiT and sinhtu^T in (25), which we write as in 

Ml7fi!i = ^,[^y\n, + lY\ Mil\TJ = ^in^+A + l), (30) 

where A is angular momentum projection on the x axis. Note that v^l2 = 
3.46 ~ TT and transverse spectrum is very close to the flux tube excitations, 
while the longitudinal found in 6J is new. 

The resulting spectrum ()3Up is in a good agreement with lattice calcula- 
tions 0, see also discussion in 0. 
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b) We now turn to the case b) aR^ <^ 1, and from (j26|) find that v 
fQ = \-^) and the equivalent of Eq. (j^ is 



G{x, y) = Dfi exp 



/i 3(3(t)2/3 ^ aR^ fcxfi\y^' 



2 2 /iV3 2 V 9 



T . (31) 



To check accuracy of our approach we can calculate the hybrid mass in the 
limit R —>■ 0, which coincides with the gluelump case [TJ. Defining for R -^ 
the stationary point jj, = Hq from the exponent of (jHTj) one has jjl = Hq = v^3a 
and the gluelump mass is 

J-^J- gluelump ^ 2,y OG \'-'^) 

which should be compared with the dedicated gluelump calculations in [7]: 

Mq = 2 f Ij (2crarfj) = 2-\/3.096(T, where we have used the value of the 
first zero of the Airy function, a = 2.338. One can see agreement on the level 
of 1.5%. 

Taking into account the last term in the exponent of ()31|) one obtains the 
lowest hybrid mass at small aR^ <C 1 



aR 



Mhybrrd{R) = 2V3a + ^"^3 + 0(i?^) (33) 

which coincides with the mass spectrum obtained in [H] for this limiting case 
by a different method. 

Thus our approach can be used as a good zeroth order approximation for 
the confined gluon Green's function and its spectrum. 

In what follows we shall use the dependence /io(-R) = /i given above in 

two limiting cases: 

a\i/3 



/io(i?) = ( -^ ) , (tR' » 1 (34) 



/io(i?) = V3a, aR' < 1 

5 The confined Coulomb interaction 

In this section we consider the confined gluon Green's function for the initial 
and final conditions corresponding to the Coulomb gluon exchange. 



With the same notations as in (|22j ) - (j24j) . one has 

GiR, T)= I DvDv^e-^Gf^ {R, T, u, u, /i) (35) 

where T is the same as in (J23|l . but now G^ is not given by p5|l . but has 
another form, due to different initial and final condition: x = (0, 0, 0, 0), y = 
{R, 0, 0, T) and the same simple general formula of Appendix 1, yields 



G 



(c) 

3 



fiUJi \ 1/2 {fi + Ji + J2)UJ± 



27rsinhc(JiT/ 27rsinh(ct;_LT) 



expJ ^^:^p——\R^coshuJiT + 2R^^^(l-coshiUiT)]\ (36) 

[ 2 smh(ci;iT) u + u J 

where we have defined u± as in previous section, while 



V/i(z/ + z/) ^„^ J27rn{n) 



so that / D/i exp I — I /q /i(t)(it|- = 1. Minimizing in (z/ — z/), one obtains 
z/ = z/, wi ^ u;f ) = .If. uj± -^ uf^ = J ,,,llr. and one has 






where 



To = -T + ^ + i lnsinh(c.rT) - ^ ln(^c.r ) + '^^{u.'?T) + 

+ In sinh(wf T) - In ((/i + Ji + J2)4°0 ' (38) 

and 

x(l + coshx) X 

"^^^^ ^ — 2sinhx — ' '^^^^ " ' ^^"^ ^ °°) ~ 2 ^^^' 

To proceed one should find the stationary point of Fq with respect to 






= 0. This is easy to do at large i?, since then z/q ~ # and 

v=vo ^ 



R 

(0) ^ [^ ^ '-'^"^ 
do not contribute. 



u;} = J^ ^ 0, cjjl^ ^ so that the last three terms on the r.h.s. of ()38p 
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and 

,....(f.)-.(..f,^). («) 

To get the modified Coulomb interaction at large R, one considers the 
integral 

/oo /"oo 9 poo 

dTG{R,T) = 2 dTG{R,T)= / e-^^ciT. (42) 

-oo Jo [iTiyi'^Zll Jo 

Inserting ()41|) into ()42|) one obtains 



y*(i?) 



27r3/2 



sinh (^) [sinh (^T) 



1/2 



Introducing new variable r = ^|r this integral can be reduced to the form 

y.(r,. v~^{m' + ^) f-dr e-K^) _m 

^ ^ 2.3/2^ J^ -^(sinhA)^/%inh(^)"4-i? ^ ^ 



where we have defined A = {a^iR^Y^'^—^B.^ooi.crR^y/'^. Finally ^{R) can be 
written as 

am = \/f 2A (^1 + ^) /(A), A = (ai^2)2/^ (45) 

/(A) = / c^y ^ (46) 



sinh (Ay) sinh(A/3A?/) 

For A ^ one has /(A) ~ ^v f? and ,^(A -^ 0) is close to unity. The exphcit 
behaviour of ^(-R) is given in the Table. 
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Table : The screening factor C,{R) and /(A) as functions of distance 
R 



A 


0.1 


0.2 


0.4 


0.6 


0.8 


1.0 


2.7 


5.4 


/(A) 


4.60 


1.91 


0.74 


0.41 


0.26 


0.18 


0.024 


0.0029 


R{im) 


0.084 


0.14 


0.23 


0.32 


0.39 


0.46 


0.98 


1.645 


^R) 


0.929 


0.796 


0.656 


0.577 


0.515 


0.469 


0.241 


0.086 



1/3 



For aR^ —>■ oo, one has from (jH^ f^ ~ iji) ^^'^ inserting this into (44) 
one obtains the asymptotics. 



V*{R) 




2 Aexp(-|) 

f^ 

R 



V2{aR 



■2\2/3 



exp 



-UaR 



2\2/3 



Vsn^/m 



(aR^) > 1 



(47) 

One can see in the Table that the screening starts at rather small values of 
R and at i? ~ 0.5 fm, the coefficient C,{R) ~ 0.5. 



6 Conclusions 

The overall static interaction between Q and Q can be derived from Eq. (jH)), 
where the term {W^^\B)) gives the confining term Kon/(-R); while the second 
term on the r.h.s. of (0) provides the screened Coulomb potential Vq{R), so 
that one has for i? ^ A, A ~ 0.2 fm. 



VstaUciR) = (tR+VS{R),VS{R) 



02^3 

R 



^R), 



(48) 



where ^{R) is given in (45), (46) and in the Table. 

In (P5|) the interference of perturbative field a^ and nonperturbative 5^ 
is taken into account. At small R,R < A, there exists another interference 
effect, which was treated before in ^2], and which provides linear behaviour of 
Vconf{R) at very small R, while without this interference Vconf{R) ~ constR^, 
for i?< A [in]. 

The behaviour aR — -| was checked on the lattice in the interval 0.1 fm 
< i? <1 fm with good accuracy Jl] and the region 0.8 fm < i? <1.5 fm 
was also measured ^^1 in the regime where string breaking is expected. Re- 
cently a thorough analysis of the region i? < 1 fm [16] has revealed that 
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^{R) is approximately constant in this interval and coincides within 15% 
with the bosonic string Casimir energy prediction (see [16] for discussion and 
earlier references). Also the heavy quarkonia spectrum calculated assum- 
ing ^(-R) = 1 is in good agreement with experiment including high excited 
bottomonium states [17,18] and deviation of ^(-R) from unity for R > 0.5 
fm seems to deteriorate this agreement. Thus the calculated in this paper 
screening factor ^{R), Eq.(45) and Table, is in a serious conflict with lattice 
and phenomenological (experimental) data. A possible solution of this para- 
dox lies in adopting the bosonic string term [16] at distances beyond 0.5 fm, 
so that the sum of the screened Coulomb term and the bosonic string term 
could imitate the original unscreened Coulomb interaction. This picture of 
trasmutation of the Coulomb into the string vibration term, if realistic, can 
be supported by the Casimir scaling study of the Coulomb-like term at dis- 
tances around 1 fm. The accuracy of the previous study by Bali [14] was 
insufficient to draw definite conclusions about the presence of the bosonic 
string term in this region. 

In addition one should recalculate the bosonic string term using the re- 
alistic hybrid spectrum found in [6], which will be reported elsewhere. 

The author is grateful for useful discussions to A.M.Badalian, N.O.Agasyan, 
A.B.Kaidalov, Yu.S.Kalashnikova and V.I.Shevchenko. 

The work is supported by the Federal Program of the Russian Ministry 
of industry. Science and Technology No. 40. 052. 1.1. 1112. 



Appendix 1 

Derivation of the gluon Green's function in the toy model of the 
confining plane (chapter 3) 

Eq.(19) describes Green function of free motion in the plane (xijX^) and 
factorizable motion in oscillator potential in the directions X2 and 2:3. For 
the latter one can use the standard textbook formula (see e.g. [121) for the 
Green's function G{xa,ta;Xb,th) corresponding to the Lagrangian 
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which we write in the Euchdean metrics 



mujQ \ ^/^ 



27rsinhc(JoT, 

exp {~ 2sinh(!!> T ^*^^" ^ '^^'* cosho^oT - 2xaXb]j 

where T = ta — h. Changing mujQ -^ ^, m = |, and identifying Xa = Xh = 0, 
one arrives at the result written in Eq. (20). 

Appendix 2 

The gluon propagator in the einbein path-integral representa- 
tion 

One starts with FSR for the free propagator which can be written as 

("OO 

G{x,y)= dsiDz)^yexp{-K) (A2.1) 

and introduce the einbein variable-dynamical mass /i(t) as in [10], so that K 
can be rewritten as 

4io V dr J 

"rft/^ + ^ + ^f^l }. (A2.2) 

[ 2^(t) 2 2 \ dt ' ' ^ ' 

In {Dz)xy, there is integration over time components of the path, namely 

('D^4)-n^|^*(i:A.4-r) (A2.3) 

where T = x^ — y^, and using the definition of fj.{t),ii{t) = ^^ one can 
rewrite the integration element in ()A2.3|) as follows {t = z^) 



dAzJn) , , . fs du(n)\'At r- /At /.„ .n 

^^ ' 2dfi{n)J-= ^y , v^=J-— -. (A2.4) 



V4^ VTT ^27ifi{n) V2/i(n) 
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Moreover the 5-function acquires the form 

6 (J2 A-24 -T) = S{s2fi - T) (A2.5) 

where we have defined 

/i= - r 2^i{T)dT. (A2.6) 

S JO 

As a resuh one can integrate in ()A2.1jl over ds using 5- function ()A2.5j) . and 
rewrite ds{Dz)xy as ds{D'^z)xy = {D^z)x.yDfj,. One can write the Green's 
function as follows 

where K is given in ()A2.2|1 . and /, l^ are 

The integration over (i'^Azj(r2) yields 

'^f^-^' = /|^^-""-^'-^"'*""('"'^)^(fl^). (A2.8) 

Taking into account the integral 

r^ di4nl^-^[,in)^Pl±f) ^ /2^^-A*v/p^W .A2 9) 

^0 y^ V At 

one has finally 



where we have used relation following from the stationary point in the integral 

fi= r fi{T)dT = ^p2 + m2. (A2.11) 

This can be compared with the integral 

which reduces to flA2.10|) after integrating over dp^ for T > 0. Eq. ()A2.12|1 
is the standard form of the free propagator. Now in the case of interacting 
gluon as in section 4 one can still use ()A2.5|) . ()A2.6|) with the result, that 
/i = /io, with fiQ defined in (34) since fiQ does not depend on time. 

15 



References 

[1] B.S.De Witt, Phys. Rev. 162, 1195, 1239 (1967) 
J.Honerkamp, Nucl. Phys. B 48, 269 (1972); 

G.'t Hooft Nucl. Phys. B 62, 444 (1973), Lectures at Karpacz, in : Acta 
Univ. Wratislaviensis 368, 345 (1976); 
L.F.Abbot, Nucl. Phys. B 185, 189 (1981). 

[2] Yu.A.Simonov, Phys. At. Nucl. 58, 107 (1995); Yad. Fiz. 58, 113 (1995). 

[3] Yu.A.Simonov, in: Lecture Notes in Physics, v. 479, P-144j Springer, 
1996. 

[4] A.B.Kaidalov and Yu.A.Simonov, Phys. Lett. B 477 163 (2000); 
hep-ph/991 2434, 

A.B.Kaidalov and Yu.A.S imonov, Phys. A tom. Nucl. 63, 1428 (2000); 
Yad. Fiz. 63 1428 (2000); |hep^h/9911291| 

[5] Yu.A.Simonov, Nucl. Phys. B (proc. Suppl) 23, 283 (1991); 

Yu.S.Kalashnikova, Yu.B.Yufryakov, Phys. Lett.B 359, 175 (1995). 

[6] Yu.S.Kalashnikova, D.S.Kuzmenko, Phys. Atom. Nucl. 64, 1716 (2001); 
ibid. 66, 955 (2003), hep-ph/0203128, hep-ph/03022073 . 

[7] Yu.A.Simonov, Nucl. Phys. B 592, 350 (2001); |hi^h/0003114[ 

[8] A.Di Giacomo and H.Panagopoulos, Phys. Lett. B 285, 133 (1992); 
A.Di Giacomo, E.Meggiolaro and H.Panagopoulos, Nucl. Phys. B 483, 
371 (1997); 

M.D'Elia, A.Di Giacomo and E.Meggiolaro Phys.Lett. B 408, 315 
(1997). 

[9] K.J.Juge, J.Kuti and C.Morningstar, Nucl. Phys. Proc. Suppl. 73, 360 
(1999); hep-lat/9809015, Phys. Rev. Lett. 90, 161601 (2003). 

[10] A.Yu.Dubin, A.B.Kaidalov, and Yu.A.Simonov, Phys. Lett. B 323, 41 
(1994); 
A.Yu.Dubin, A.B.Kaidalov, Yu.A.Simonov, Phys. At. Nucl. 56, 1745 



(1993); plep^ph/9311344 



16 



[11] Yu.A.Simonov and J.A.Tjon, Ann Phys. 228, 1 (1993); 

Yu.A.Simonov, J.A. Tjon, Annals Phys. 300, 54 (2002); 
|hi3-ph/0205165| 



[12] Yu.A.Simonov, Phys. Kept. 320, 265 (1999); JEPT Lett 69, 505 (1999). 

[13] H.G.Dosch, Phys. Lett. B 190, 177 (1987); 

H.G.Dosch and Yu.A.Simonov, Phys. Lett. B 205, 339 (1988); 
Yu.A.Simonov, Nucl. Phys. B 307, 512 (1988). 

[14] G.Bah, Nucl. Phys. Proc. Suppl. 83, 422 (2000), Phys. Rev. D62, 114503 
(2000); 

A.Hasenfratz, R.Hoffmann and F.Knechtli, Nucl. Phys. Proc. Suppl. 
106, 418 (2002). 

[15] B. Bolder, T.Struckmann, G.S.Bali, N.Eicker et al., Phys. Rev. D63, 
074504 (2001). 

[16] M.Liischer and P.Weisz, JHEP 0207, 049 (2002). 

[17] E.Eichten et al. Phys. Rev. D 21, 203 (1980); 

E.J.Eichten, K.Lane and C.Quigg, Phys. Rev. Lett. 89, 162002 (2002). 

[18] A.M.Badalian, B.L.G.Bakker, Phys. Rev. D 67, 071901 (2003); 

A.M.Badalian, B.L.G.Bakker and A.I.Veselov, Phys. Rev. D (in press), 
[hep-ph/0311010. 

[19] R.P.Feynmann, A.R.Hibbs, Quantum Mechanics and Path Integrals 
McGraw-Hill, N.Y., 1965, problem 3.8. 



17 



